Mesoscopic effects in superconductor-ferromagnet-superconductor junctions 



A.Yu. Zyuzin \ B. Spivak 2 , M. Hruska 2 

1 A.F.Ioffe Physical- Technical institute, 194021 St. Petersburg, Russia 

2 Physics Department, University of Washington, Seattle, WA 98195 

We show that at zero temperature the supercurrent through the superconductor - ferromagnetic 
metal - superconductor junctions does not decay exponentially with the thickness L of the junction. 
At large L it has a random sample-specific sign which can change with a change in temperature. In 
the case of mesoscopic junctions the phase of the order parameter in the ground state is a random 
sample-specific quantity. In the case of junctions of large area the ground state phase difference is 
±tt/2. 



This work has been motivated by recent experiments on the superconductor-metallic ferromagnet- 

superconductor junctions where the ground state of the system with the superconducting phase difference equal 
to 7r has been observed. It is well known that the sign of the critical supercurrent of pure SFS junctions oscillates 
with the width L of the ferromagnetic region |^-|^]. This is due to the difference in Fermi momentum of electrons 
with different spins at the Fermi energy causing the superconducting wave function in the ferromagnetic region to 
oscillate with the characteristic distance h/(k^ — k\). Here fcj, k± are Fermi momenta of up and down spins, which 
are different because of the finite exchange spin-splitting energy / in ferromagnets. 

It is important to note that in pure junctions the characteristic distance of oscillations is inversely proportional to 
I and that at zero temperature the modulus of the critical current does not decay exponentially with L. 

In the opposite limit of disordered ferromagnets L > Lj = ^J^- > I the average critical current (7 C ) decays 

exponentially with L and with temperature T (lj] 



(where brackets () stand for averaging over the random realizations of the scattering potential and D is the the 
electron diffusion coefficient in the ferromagnet. It is interesting however, that in addition to the exponential decay, 
the average critical current oscillates as a function of L. It also changes sign as a function of T, provided D/L 2 < I. 

The value of (I c ) can be negative which means that in the ground state of the junction the phase difference <f> of 
the superconducting order parameter is tt rather than zero. The experimentally measured critical current in this case 
is the absolute value |7 C |. 

The oscillations of the average critical current have been observed experimentally |f|,||. We would like to note, 
however, that according to Eq.l , these oscillations can be observed only in the case when the exchange spin-splitting 
energy in the ferromagnet is relatively small so that the characteristic distance of oscillations is large. This limits 
significantly the choice of ferromagnetic metals which can be used in the junctions to observe this effect. 

In this paper we show that the exponential decay of the supercurrent with L in Eq. 1 originates from the averaging 
procedure. Before averaging over the impurity configurations at L Lj the supercurrent I s (cj>,L) has a random 
sample-specific sign while its modulus does not decay exponentially. The practical consequence of this conclusion is 
that the Josephson effect survives in the case of ferromagnets with large I when L ^> Lj. We would like to mention 
that this feature is a particular case of a general statement that the Friedel oscillations in disordered metals do not 
decay exponentially at zero temperature . 

Below we discuss the limit of thick SFS junctions L Li when their superconducting properties are determined 
by mesoscopic effects. We show that in this case the critical current of the junctions does not depend on the spin 
splitting /. The ground state of such a junction is, generally speaking, doubly degenerate with the phase difference 
having a random sample-specific modulus distributed in the interval (0,7r). We also show that the critical current 
undergoes random oscillations as a function of T. Thus the junctions with L Lj should exhibit the same sequence 
of effects as the junctions with L < Lj. 

The energy of the Josephson junction E{<f>) is an even and periodic function of the phase of the order parameter. 
It can be represented in the form 

oo 

E{(j)) = J2 E n cos n<f> (2) 

n=l 



f 



while the current through the contact is determined by the relation J = 2e^. The coefficients E n are random sample- 
specific quantities. For L> Li all average coefficients (E n ) ~ exp(— L/Lj) are exponentially small. In this case the 
typical values of the coefficients E n can be estimated from their variances. To simplify the analysis we consider the 
case when the superconductors and the fcrromagnet are separated by tunneling barriers of small transparency. We 
will show below that at D/L 2 >T and L > £ we have 



(where S is the area of the junction, is the density of states per spin in the ferromagnet, £ is the zero-temperature 
coherence length in the superconductor and g is the conductance per unit area of the surface between the ferromagnet 
and the superconductor. Eqs. 3,4 correspond to the diagrams shown in Figures l.b),c). In the case when the 
tunneling transmission coefficient of the insulator between the superconductor and the ferromagnet is small, the 
typical amplitude (((E 2 ) 2 )y/ 2 of the second harmonic contains an additional power of g compared to the amplitude 
(((Ei) 2 )) 1 / 2 of the first harmonic and therefore for typical samples of low g, E(<f>) is well approximated by the first 
harmonic. Since (((Ei) 2 )) 1 / 2 does not decay exponentially with L, it is much larger than (Ei). This means that E\ 
has a random sign. In the cases when E\ < the ground state of the junction corresponds to <Pqs = it. 

The correlation function (E\E 2 ) contains an additional power of g with respect to ((E 2 ) 2 ) and can be neglected 
which means that Ex and E 2 are random uncorrelated quantities. In rare samples where the amplitude of the first 
harmonic is small (E\ <~ E 2 ) one has to take into account the second harmonic in Eq.2. In this case the ground state 
is doubly degenerate and corresponds to the phase difference 

<j) GS ~ ± arccos (5) 

4£/2 

Therefore the sample-specific random absolute value of the ground state phase difference is in the interval < \<j>as\ < 
ir/2, provided |£i| < 4|_E 2 |- The probability of such an event is of order 



'((E2) 2 ) 



£L(-^-) 2 (6) 



\/ ((£i) 2 ) y 8irv Q D 

We would like to mention that in the case of a transparent superconductor-ferromagnet boundary (j)cs is a random 
sample-specific quantity of order one. 

Let us now discuss the temperature dependence of the critical current I C (T). To do so we calculate the correlation 
function 

{E\(Ti)Ei(T 2 )) ~ exp(-2i v / 7 r(T 1 +T 2 )/Z?) (7) 



It follows from Eq.7 that the quantities ((E^T)) 2 ) ~ cl^-p(-^ /2 L^/T/D) and ((Ei(T)Ei(O)) ~ cxp(- 2L^/TjD)) 
decay with different rates as T increases. This indicates that in addition to an exponential decay, the quantity Ei(T) 
exhibits random sample-specific oscillations with a period of order T* ~ D/L 2 . 

The results presented above were obtained in the approximation when the variations of the phase of the order 
parameter along the superconductor-ferromagnet surface are neglected. This is a good approximation for the samples 
of small area. Below we show that in the samples of large area the possibility of spatial fluctuations of the order 
parameter phase along the superconductor- ferromagnet surface leads to an average critical current which is propor- 
tional to the area of the junction and does not decay exponentially even in the limit when L 3> Li. In this case the 
ground state of the junction is doubly degenerate and 

(Pas = ±tt/2 (8) 
This can be shown using an expression for the effective Josephson energy 



= J E c(p) cos 4>(p)dp 



E = / E c (p) cos (j)(p)d P + — 1 / dr(V0(rr) (9) 
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where p is the coordinate along the surface between the superconductor and the ferromagnet, the r-integration is 
performed in the bulk of the superconductors, m is the electron mass, and N s is the superfluid density in super- 
conductors. Eq.9 is valid on a scale larger than L along the surface. The random function E c (y) is character- 
ized by its average (E c (p)) = and a quickly decaying (at \p — p'\ > L) correlation function (E c (p)E c (p')) ~ 
{E c (p - p')E c (0)) ~ ((Ei(S = L 2 )) 2 )/L 4 . The phase difference <j>{p) = (0) + 5<f>{p) is a random function of p. At 
small Ei 6cj)(p) <C 1. Minimizing Eq.9 with respect to S<p(p) we get an expression for the effective energy per unit 
area || 

E eff = -E e 1 ff sin 2 {cl>) (10) 

(where E\^ ~ TTT[l? ((Ei(S = L 2 )) 2 )), determining the average critical current density as (I c ) = 2eE\^ and the 
average phase difference in the ground state as given by Eq.8. 

The results in Eqs.3,4 were calculated microscopically describing the SFS junction by a Hamiltonian of the form 

H = Hbcs + Ht + Hp (11) 

where Hbcs is the Hamiltonian of superconducting leads, the Hamiltonian 

R T = t j dp i (*f(p i ,Zi;a)¥ F (Pi,Zi;a) + h.c) (12) 

i=l, 2:q 

describes tunneling between superconductors labeled by i = 1, 2 and the ferromagnetic metal {z\ = 0, z% = L), labeled 
by index F, and a is the spin index. The integration is taken over the surfaces between superconductors and the 
ferromagnetic metal. The last term in Eq.ll corresponds to the disordered ferromagnetic metal unperturbed by the 
presence of superconductors, where spin up and down bands are split by the exchange field / 



H F = H Q + I J drf F (r; a) a z a ^ F (r; (3) . (13) 

Here Hq is the Hamiltonian of noninteracting electrons which contains the operators of the kinetic energy and a 
random field U (r) . We assume that the random potential is white-noise correlated so (U (r)) = and (U (r) U (r')) = 
{2itvqt) S (r — r') (where r is the mean free time). 

In the lowest order in tunneling through the superconductor-ferromagnet boundary we get 

Ei = — V t 4 [ dp 1 dp 2 dp' 1 dp' 2 {F + (e fe ; p x , 0; pi, 0) G_ Q; _ Q (e fe ; p[, 0; p' 2 , L) F (e fe ; p' 2 , L; p 2 ,L) G a - a (e fc ; p 2 , L; p 1 ,0) + h.c} 



c k :a 



's i; s 2 

(14) 

where = n(2k + 1)T is the Matzubara frequency and k = 0, 1, ... is an integer. We use the usual diagrammatic 
technique for averaging the products of electron Green's functions pOj] . Diagrams for correlation functions ((£i) 



and (^(E 2 ) 2 ^ are shown in Fig.l.b),c). There are two important blocks in these diagrams. The first one corresponds 

to diffuson and cooperon ladder made of single particle Green's functions in the ferromagnet. In the case of a large / 
only diffusons and cooperons for parallel spins survive. They are equal to each other and the same as in the absence 
of the external magnetic field, satisfying the equation |ll[ 

(-DV 2 + \e k - e k ,\) C{\e k - e k ,\ ; r, r') = 9 {-e k e w )8 (r - r') (15) 

We would also like to mention that for It > 1 we can neglect contributions of diagrams shown in Fig.l.d). 

The second block is the ladder made of the anomalous Green's functions F(e k ;r,r') in the superconductor. This 
average can be expressed in terms of the averaged product of the advanced and retarded Green's function in the 
normal metal. Since 

oo 

F(e k ;v,v') = -i J (G R (E, r, r') - G A (E, r, r')) (16) 

— oo 

(where fl 2 = A 2 + e 2 ,), the averaging (F(e k ; r, r')F(e' k ; r, r')) is equivalent to the averaging of Green's functions in the 
normal metal. The diffusion propagator (F(efc;r,r')F(e^.;r,r')) ~ P(Sl,il';r,r') obeys the diffusion equation 
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(~d\7 2 + n + n>) p(n, n' ; r, r') = i£ A* (r _ r ') (17) 

The result of integration of four electron Green's functions over the surface in diagrams shown in Fig.l.b),c) is 
estimated as r 2 g , where g is the dimensional tunneling conductance per unit area of the boundary. 

The solution of Eq.15 satisfying the boundary conditions — Cflefe — e^/l ;r, r') = at the superconductor - ferro- 

dz 

magnet metal surfaces is given by ( for z > z' ) 

c(| £fc - e ,|;ry) = *(-^) / cosh(Q(z - L ^l^ {z ' + £/2) (18 ) 



where f / = y 1/ 2 + ■ fc — fc ■ . Substituting this expression into diagrams shown in Fig.l.b)c) we get 



(ft) 2 



S ^^^.?. e (--0/|^^-i W (20) 



Calculating integrals in Eqs. 19,20 we arrive at Eqs.3-4. 

To calculate the temperature dependence (E\ (T\) E\ (T 2 )) we should substitute = 7r(2fc + l)Ti and = 
Tr(2k + 1)T 2 into Eq.19 . Then for L Tl = (D /T^) 1 ' 2 ,L T2 = {D/T^ 1 / 2 < L we get (sinhQL) -2 = 4cxp[-2L(q 2 + 
(ttTi + 7rT 2 )/L>) 1/2 ] and finally we arrive to Eq.7. 

In conclusion we have shown that the critical current of a mesoscopic SFS junction at small temperatures does not 
decay exponentially with the ferromagnet thickness. It has a random sign, which changes with temperature. The 
ground state phase difference of the junction is a random quantity < (j)cs < 7r/2. In the case of junctions of large 
area the phase difference is c/)gs = tt/2. Let us estimate the typical value of the critical current using Eq.3. Taking, 
for example, iron as a ferromagnet with L ~ £ = 3 x 10~ 6 cm, the area of the surface S ~ 10~ 7 cm 2 , the elastic 
mean free path I ~ 10~ 6 cm and assuming that the transmission coefficient through the superconductor-ferromagnet 
boundary is of order one we get I c ~ 10 -6 A. The estimate based on Eq.3 scales with the junction area as y/S and is 
valid at small 5*. For junctions of large area, Eq.10 implies that the critical current is proportional to S. 
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FIG. 1. The diagrams contributing to: a) the current, to lowest order in transparency, b) the correlator (EiEi), c) 
the correlator (E2E2), d) a diagram insignificantly contributing to the correlator (E\Ei) in a superconductor - ferromagnet 
- superconductor junction. The symbols Si and S2 correspond to the first and second superconductor respectively and F 
denotes the ferromagnet. Thin lines represent the single electron Green's functions, thick lines represent the anomalous Green's 
functions, dashed lines denote averaging over the configurations of the impurity scattering potential. Ladders formed of Green's 
functions and impurity scatterings in the ferromagnet correspond to diffuson or cooperon ladders. 
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